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^ SUMMARY 
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the  contact  pressure  distribution  is  then  determined  by  an  iterative  procedure  which  allows 
for  compatibility  of  displacement  between  pin  and  lug.  Having  found  the  contact  pressure 
the  stress  distribution  throughout  the  lug  is  determined  by  a routine  finite  element  analysis. 
This  method  is  applicable  to  cases  of  neat,  clearance  or  interference  fit  pins.  Numerical 
results  for  neat  fit  pins  show  reasonable  agreement  with  some  test  data. 
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NOTATION 

{S/(}  Vector  of  radial  displacements  at  boundary  nodes 

(R  Radial  strain 

fT  Tangential  strain 

OR  Radial  stress 

OT  Tangential  stress 

ORT  Shear  stress 

V Poisson’s  ratio 

0 Angular  co-ordinate 

[A]  Matrix  of  Influence  Coefficients  relating  radial  loads  to  radial  displacements 

{Pfl}  A vector  of  radial  loads 

{P}  Initial  load  distribution 

{Pm}  Modified  load  distribution 

[XA]  Matrix  relating  X displacements  of  nodes  to  radial  loads 

lYA]  Matrix  relating  T displacements  of  nodes  to  radial  loads 

{XCOORD}  Matrix  of  X Co-ordinates  of  node  points 
{YCOORD}  Matrix  of  Y Co-ordinates  of  node  points 
DISPXi 
DISPYi 
E 
F 

ICL 
K 
P 
R 

RDISii 
RH 
RLUG 
RPIN 
d 


J 


Displacement  in  X direction  at  node  i due  to  load  distribution 
Displacement  in  Y direction  at  node  i due  to  load  distribution 
Young’s  modulus 
Fit  of  pin  in  lug 

Node  number  defining  extent  of  contact  arc 
Stress  concentration  factor 
Pin  load 
Radius 

Radial  displacement  of  node  i due  to  a unit  radial  load  at  node  j 

Radial  vector  to  displaced  node  on  lug  hole  boundary 

Radius  of  Lug 

Radius  of  pin 

Vertical  translation  of  lug 

Thickness  of  lug. 


1.  INTRODUCTION 


The  pin-loaded  lug  shown  in  Figure  1 is  a common  form  of  connection  in  aircraft  structures. 
Until  recently,  the  stress  analysis  of  a lug  was  usually  carried  out  in  a rather  simple  fashion 
consideration  usually  being  given  only  to  its  ultimate  strength  which  was  assessed  by  assuming 
simple  stress  distributions  possibly  modified  by  empirical  factors.  However,  when  considering 
the  fatigue  performance  of  a lug,  or  its  strength  in  a cracked  condition,  an  accurate  knowledge 
of  the  stress  distribution  is  required. 

When  the  contact  pressures  applied  by  the  pin  to  the  lug  are  known  then  the  stress  analysis 
of  the  lug  can  be  carried  out  using  the  finite  element  method.  It  has  been  often  assumed  that 
the  contact  pressure  p,  for  a lug  with  the  force  applied  along  the  centre  line,  is  given  by 

p = (2/’/7c/?f)  cos  0 1^1^  (I) 

p = 0 I ^ I > 

where  P is  the  Pin  load 

R is  the  radius  of  the  pin 

t is  the  thickness  of  the  lug 

0 is  the  angular  Co-ordinate  shown  in  Figure  2. 

Intuitively,  it  might  be  expected  that  the  actual  contact  pressure  would  resemble  that  given 
by  equation  (1).  However  equation  (I)  has  two  important  defects: 

1)  It  does  not  include  the  influence  of  flexibility  of  the  pin  and  lug  on  the  contact  pressure. 

2)  The  area  of  contact  between  pin  and  lug  is  assumed  fixed,  no  allowance  being  made 
for  changing  areas  of  contact  between  pin  and  lug  as  a result  of  the  type  of  fit  or  size 
of  load  applied  to  the  pin. 

Also,  experimental  results,  reference  (I),  have  shown  that  the  lug  stresses  and  strains  vary 
non-linearly  with  pin  load  whilst  equation  (I)  predicts  a linear  variation. 

A method  for  determining  the  contact  pressures  for  the  case  of  flexible  lug  loaded  by  a 
rigid  pin  is  described  in  references  (2)  and  (3).  A procedure  in  which  the  flexibility  of  both  the  lug 
and  the  pin  is  considered  is  given  in  reference  (1).  Both  of  these  methods  use  the  finite  element 
approach.  Reference  (4)  gives  details  of  a photo-elastic  investigation  and  provides  some  data 
which  can  be  compared  with  theoretical  results. 

The  present  work  arose  out  of  the  need  to  develop  an  accurate  method  for  lug  stress  analysis 
which  could  be  used  in  conjunction  with  an  existing  general  purpose  finite  element  program 
described  in  reference  (5).  The  method  is  similar  in  principle  to  that  described  in  references  (2) 
and  (3)  but  differs  in  the  details  of  the  computer  implementation.  It  is  also  restricted  to  the  case 
of  a flexible  lug  loaded  by  a rigid  pin.  This  is  considered  a reasonable  approximation  for  an 
aluminium  alloy  lug  loaded  by  a steel  pin  which  is  a common  situation  in  aircraft  structures. 
Also,  it  is  assumed  that  frictional  forces  between  the  pin  and  the  lug  can  be  ignored  and  that 
the  problem  can  be  treated  as  a two-dimensional  (plane  stress)  one.  A particular  lug  is  considered 
here  for  the  case  of  neat,  clearance  and  interference  fit  pins.  For  the  neat  fit  pin  the  theoretical 
results  are  compared  with  some  strain  gauge  results. 

2.  METHOD 

2.1  Description  of  the  lug 

The  method  is  described  by  its  application  to  the  particular  lug  shown  in  Figure  3;  this  is 
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of  7079-T6  aluminium  alloy  and  formed  part  of  a rudder  actuator  fitting  of  a large  transport 
aircraft.  Relevant  properties  of  the  alloy  lug  were  taken  to  be 

Young’s  modulus  E = 71-02  x 10’ MPa 

Poisson’s  ratio  v = 0-32 

The  line  of  action  of  the  pin  load  is  taken  to  lie  along  the  centre  line  of  the  lug.  The  straight 
end  of  the  lug  is  constrained  in  a way  that  leads  to  a uniform  stress  being  applied  at  that  end. 
Since  the  structure  is  symmetric  about  the  centre  line  only  half  of  it  need  be  considered  in  the 
^ analysis.  The  finite  element  idealization  of  the  lug  based  on  linear  strain  triangles  is  shown  in 

^ Figure  4;  also  shown  is  the  cartesian  co-ordinate  system  which  has  its  origin  set  at  the  centre 

of  the  lug  hole.  On  the  inside  of  the  lug  hole  are  a series  of  boundary  nodes  at  which  forces 
and  displacements  between  the  pin  and  lug  are  to  be  considered.  These  nodes  are  numbered  1 
. to  17.  The  method,  involves  an  iterative  procedure,  which  determines  the  forces  required  at 

these  boundary  nodes  to  produce  prescribed  displacements;  this  is  achieved  by  setting  up  a 
matrix  of  influence  coefficients  which  defines  the  flexibility  of  the  lug.  In  the  following,  as  far  as 
possible,  the  notation  is  chosen  to  agree  with  that  used  in  the  computer  program. 

2.2  Matrix  of  Influence  Coefficients 

I 

, The  radial  force  applied  at  boundary  node  j is  denoted  by  PriU  = 1,  NPOINT);  this  is 

taken  to  be  positive  when  applying  pressure  on  the  lug.  Also  the  x and  y displacements  at 
boundary  node  i are  denoted  by  DISPXi  and  DISPY t (i  = 1,  NPOINT).  The  following  matrix 
equations  relate  these  quantities: 

{DISPX}  = IXA]  {Pr}  (2) 

{D/SPT}  = (3) 

where  [XA]  and  [YA]  are  square  matrices  of  order  NPOINT  x NPOINT.  A typical  element, 
XAi},  of  the  first  matrix  is  simply  the  value  of  DISPXi  for  Pr]  = 1 with  all  other  values  of  Pr 
equal  to  zero;  an  analogous  interpretation  holds  for  YAij.  These  matrices,  [XA]  and  [YA], 
can  be  determined  from  a finite  element  analysis  for  N POINT  \oad  cases,  each  case  corresponding 
to  a unit  radial  load  at  one  node. 

It  is  more  convenient  to  define  a single  matrix  [A]  relating  radial  displacements  Sri  (taken 
positive  when  directed  outwards  from  the  origin)  to  the  radial  loads;  this  is  given  by 

{«„}  = (4) 

Matrix  [A]  can  be  obtained  from  [XA]  and  [YA].  Referring  to  Figure  5,  the  radial  distance 
to  node  / after  its  displacement  due  to  a unit  radial  load  at  nodey  is  given  by 

RDIS,)=  ^{(XCOORDi  + XAij)^  + (YCOORDi  + YAi,)^}  (5) 

where  XCOORDt  and  YCOORDt  are  X and  Y co-ordinates  of  point  i.  The  yth  element  of 
[A]  is  simply  the  displacement  of  node  i for  this  loading: 

Ai]  = RDISij  - RLUG  (6) 

• where  RLUG  is  the  radius  of  the  hole. 

It  is  necessary  to  obtain  the  radial  loads  around  the  lug  in  terms  of  the  radial  displacements. 
This  requires  inversion  of  matrix  [A]  and  leads  to  the  equation 

{Pr}  = [-4]-^  {Sr}  (7) 

2.3  Initial  Load  distribution 

To  begin,  it  is  necessary  to  assume  an  initial  distribution  of  contact  pressure.  The  distribution 
given  by  equation  (1)  is  assumed  and  is  referred  to  as  the  ‘Sinusoidal  distribution’.  Replacing 
the  distributed  loads  by  discrete  radial  loads  Pt  at  the  appropriate  boundary  nodes  the  corres- 
ponding cartesian  displacements  DISPXi  and  DISPYi  at  all  boundary  nodes  are  found  from 
equations  (2)  and  (3). 
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2.4  Establishing  compatibility  in  displacement  of  Pin  and  Lug 

Generally,  the  displacements  calculated  above  will  not  be  compatible,  either  because  they 
involve  some  parts  of  the  lug  lying  inside  the  pin  boundary  or  because  they  give  a separation 
between  pin  and  lug  at  points  where  a non-zero  contact  pressure  has  been  assumed.  This  in- 
compatibility is  demonstrated  in  Figure  6,  where  the  position  of  the  lug  hole  prior  to  loading 
and  its  (exaggerated)  shape  subsequent  to  an  applied  load  are  both  shown.  Here  the  pin  is  shown 
as  a dotted  outline  and  corresponds  to  the  case  of  an  interference  fit.  The  fit  is  denoted  by  F 
and  is  based  on  the  difference  in  radii  between  pin  and  hole;  it  is  positive  for  clearance.  Hence 
the  pin  radius  is  given  by 

RPIN  = RLUG-F  (8) 

After  loading,  the  elongated  shape  of  the  lug  hole  has  undergone  a maximum  Y displacement 
DISPYi  at  boundary  node  1 which  is  at  the  vertex  of  the  lug.  This  point  is  now  at  a radial 
distance,  d,  from  the  corresponding  point  on  the  pin,  where  d is  given  by 

d = F - DISPYi  (9) 

The  deflected  shape  of  the  lug  is  translated  vertically  an  amount  d such  that  there  is  zero 
relative  displacement  between  the  bottom  of  the  pin  and  the  bottom  of  the  deflected  lug  shape. 
Consider  the  displacement  of  boundary  node  i shown  in  Figure  7 whose  co-ordinates  are 
XCOORDi,  YCOORDi.  Under  an  applied  loading,  node  / is  displaced  an  amount  DISPXt, 
DlSPYi;  after  a vertical  translation  d,  the  cartesian  co-ordinates  are  given  by: 

XCOORD,  + DISPXt  (10) 

YCOORDi  + D/SPYt  + F - DISPY i (11) 

Hence  th?  radial  distance  from  the  centre  of  the  pin  to  the  boundary  node  /'  on  the  deflected 
lug  is  given  by 

RH,  = .^{(XCOORDi  + DISPXt)^  -I-  ( YCOORDi  + DISPYi  + F - DISPY, y}  (12) 

Thus,  to  make  boundary  node  i conform  to  the  rigid  shape  of  the  pin  a radial  displacement 
Sri  must  be  applied;  this  is  given  by; 

Sr,  = RPIN-RHi  (13) 

For  the  purposes  of  determining  the  required  radial  corrections  it  is  necessary  to  establish 
the  length  over  which  the  lug  and  pin  are  in  contact.  It  is  assumed  that  the  pin  and  lug  are  in 
contact  from  node  1 up  to  a node  designated  ICL.  For  a neat  fitting  pin,  ICL  is  the  node  which 
has  angular  co-ordinates  of  tzII.  In  the  case  of  clearance  or  interference  fits,  ICL  is  the  next 
highest  node  to  that  at  which  the  current  pressure  distribution  is  down  to  zero.  Corrective  dis- 
placements given  by  equation  (13)  are  applied  to  all  nodes  up  to  ICL  and  any  higher  numbered 
nodes  for  which  is  positive.  Beyond  ICL  all  nodes  for  which  Sm  is  negative  and  thus  there 
is  a separation  between  pin  and  lug,  no  corrective  displacements  are  applied  and  Sri  is  set  to 
zero.  Note  however  that  negative  values  of  Sri  may  be  applied  within  the  contact  length. 

In  Figure  6 the  segment  AB  is  the  portion  of  the  lug  over  which  the  corrective  displacements 
are  applied.  Over  segment  BC  there  is  no  contact  between  pin  and  lug  and  hence  no  corrections 
are  applied. 


2.5  Modification  of  load  Distribution 


The  loads  that  result  from  the  corrective  displacements  are  obtained  from  equation  (7), 
and  are  added  to  the  originally  assumed  load  distribution.  No  corrective  loads  are  added  at 
points  for  which  the  corrective  displacements  are  zero.  Hence  the  modified  load  distribution 
{P\f}  is  given  by: 

{/’.w}  = [/!]-•{««}  + {/’}  (14) 

(In  the  program  the  inverse  matrix  has  been  overwritten  on  the  original  matrix.) 
Generally  this  load  distribution  will  give  a different  resultant  pin  force  to  that  required. 
Before  proceeding,  the  load  distribution  as  given  by  equation  (14)  is  scaled  to  return  the  correct 
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total  pin  force.  For  neat  and  clearance  fit  cases  the  {Pm}  are  simply  scaled  linearly  by  the  ratio 
of  the  required  pin  load  to  the  pin  load  given  by  {Pm).  This  will  not,  however,  work  for  inter- 
ference fits.  The  corrective  load  distribution  [A]  ' {8«}  for  even  small  interference  fits,  can 
produce  a large  resultant  load  in  the  opposite  direction  to  the  applied  sine  load.  If  the  magnitude 
of  the  sine  load  is  not  sufficiently  large  then  the  resultant  pin  force  will  be  in  the  opposite  direction 
to  that  required.  To  obtain  solutions  in  the  required  direction  it  is  necessary  to  scale  up  the 
initial  sine  load  until  the  resultant  pin  load  from  equation  (14)  is  just  greater  than  the  required 
pin  load.  How  close  the  resultant  load  is  to  the  required  load  depends  on  the  smallness  of 
increments  in  which  the  sine  load  is  scaled  up.  These  increments  are  usually  set  at  0 - 1 of  the 
initial  sine  load. 

This  completes  the  first  iteration. 

2.6  Further  Iterations 

The  load  distribution  just  determined  becomes  the  input  for  the  second  iteration.  The 
corresponding  deflections  are  found  from  equations  (2)  and  (3).  The  corrective  displacements 
required  are  found  from  equations  (12)  and  (13).  The  modified  load  distribution  is  found  from 
equation  (14).  After  scaling,  the  third  iteration  begins.  Iterations  are  continued  until  successive 
load  distributions  show  negligible  change. 

3.  PROGRAM 

This  program  requires  a data  file  INPUT.  The  form  of  INPUT  is  shown  in  Appendix  1. 
In  this  data  file  it  is  necessary  to  specify  the  number  of  points  around  the  lug.  the  number  of 
iterations  to  be  carried  out,  the  type  of  fit  and  the  required  pin  load  as  well  as  geometrical  data. 
In  addition,  an  initial  approximate  load  distribution  is  required.  A sinusoidal  load  distribution 
has  been  used  for  the  results  in  this  analysis;  however  a single  concentrated  load  at  the  vertex 
of  the  lug  will  give  the  same  results  (but  requiring  more  iterations)  and  is  more  convenient  to 
set  up  in  the  data.  Also  required  is  a matrix  of  influence  coefficients  which,  as  already  men- 
tioned, is  obtained  from  a standard  finite  element  analysis  of  the  lug. 

For  each  iteration,  the  radial  loads,  the  X and  Y components  of  the  radial  load,  the  X 
and  Y displacements,  the  contact  pressure  and  pin  loads  are  printed  in  file  OUTPUT.  Another 
output  file  is  CHECK;  this  is  simply  an  echo  of  the  input  data. 

A listing  of  the  program  is  given  in  Appendix  II. 

4.  RESULTS 

Distributions  of  contact  pressure  around  the  particular  lug  shown  in  Figure  3 have  been 
obtained  for  neat,  clearance  and  interference  fits  for  various  pin  loads.  Tables  I to  6 contain 
the  results  for  various  fits  for  a pin  load  of  50,000  N.  For  comparison  purposes  these  contact 
pressures  are  shown  in  Figure  8.  In  Tables  7 to  9 are  the  results  for  a clearance  fit  of  0 05  mm 
with  varying  pin  loads  from  25,000  N to  200,000  N;  these  distributions  of  contact  pressure  are 
plotted  in  Figure  9.  Results  for  an  interference  fit  of  —004  mm  for  pin  loads  varying  from 
30,000  N to  200,000  N are  shown  in  Tables  10-12  and  are  plotted  in  Figure  10. 

Returning  to  the  standard  finite  element  analysis,  the  stresses  around  the  lug  are  obtained  for 
the  cases  corresponding  to  pin  loads  of  50,000  N.  These  results  are  contained  in  Tables  13  to  18 
and  are  plotted  in  Figures  1 1 to  16.  Here,  the  stresses  are  in  the  polar  co-ordinate  system  where 
an.  OT  and  anr  denote  respectively  the  radial  and  tangential  direct  stresses  and  the  shear  stress. 

Experimental  data  in  the  form  of  strain  gauge  readings  are  available  for  this  lug  and  are 
shown  in  Table  19;  these  are  for  the  case  of  a neat  fitting  pin  with  a pin  load  of  50.(XX)  N.  The 
location  of  the  gauges  is  shown  in  Figure  17.  These  strain  gauges  have  been  placed  to  measure 
radial  strains  (n  and  tangential  strains  cr  and  correspond  to  radial  lines  through  boundary 
nodes  1 . 5 and  9 (angles  of  = 0 , 45  and  90°  respectively).  Theoretical  stresses  can  be  con- 
verted to  strains  by  the  relations: 

fH  = (<»«  - v.aT)IE 
(j.  = (ar  — v.anME 
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where  E is  Young's  modulus  and  v is  Poisson's  ratio.  These  strains  have  been  ealculated  and 
are  shown  in  Table  20.  In  Figure  18  are  shown  the  plots  for  these  theoretical  strains,  the  experi- 
mental strains  have  been  superimposed  on  these  plots  and  show  reasonably  good  agreement. 

In  reference  (4)  stress  concentration  factors  obtained  from  photo-elastic  methods  are  given 
for  various  lug  geometries.  These  stress  concentration  factors  are  for  the  tangential  stress  at 
the  edge  of  the  lug  hole  for  an  angle  of  = 90°.  In  Appendix  III,  theoretical  stress  concentration 
factors  have  been  calculated  for  clearance,  neat  and  interference  fits,  the  pin  load  being  SO.tXK)  N. 
For  the  geocr  'ry  of  the  lug  in  Figure  (3),  the  only  direct  comparison  of  results  that  can  be  made 
is  that  of  a neat  fit.  Flere  the  experimental  value  is  approximately  2 -4  and  the  trend  is  for  higher 
stress  concentration  factors  for  clearance  fits  and  lower  factors  for  interference  fits. 


DISCUSSION 

The  procedure  by  which  the  contact  load  distribution  is  obtained  in  the  program  is  an 
iterative  one.  This  is  due  to  the  non-linear  nature  of  the  contact  problem.  Solutions  tend  to 
oscillate  about  the  exact  solution  and  slowly  converge  to  it.  It  was  found  necessary  to  perform 
up  to  1000  iterations  until  two  consecutive  iterations  agreed  closely.  There  were  some  cases 
for  which  the  solution  would  not  converge  and  these  will  be  mentioned  later. 

To  examine  the  influence  of  the  type  of  fit  on  the  load  distribution  the  various  load  distri- 
butions corresponding  to  clearance,  neat  and  interference  have  been  plotted  on  Figure  8.  Firstly, 
considering  the  curve  for  the  clearance  fit,  it  is  seen  that  this  curve  exhibits  a characteristic 
peak  at  0 = 0°  which  then  drops  off  at  15°.  This  is  followed  by  a steep  increase  of  contact  pressure 
to  a maximum  at  60°;  this  then  reduces  rapidly  and  the  pin  is  free  from  the  lug  beyond  85  . 
In  the  case  of  a neat  fit,  the  peaks  in  the  curve  have  flattened  out  and  the  contact  arc  has  in- 
creased. This  trend  continues  for  the  cases  of  interference  fits.  With  the  higher  interference  fits 
the  maximum  contact  pressure  between  the  pin  and  lug  increases  as  does  the  length  of  arc  in 
contact.  For  an  interference  fit  of  — 0 07  mm  the  contact  arc  is  160°. 

Consider  now  the  effect  of  the  magnitude  of  the  pin  load  on  the  load  distribution.  As  the 
pin  load  is  increased  the  area  of  pin  and  lug  in  contact  will  change,  the  lug  tends  to  wrap  itself 
around  the  pin  to  a greater  extent.  This  is  demonstrated  in  Figure  9.  Flowever.  in  the  case  of 
an  interference  fit,  ( — 0 04  mm)  as  shown  in  Figure  10,  the  contact  length  decreases  with  increase 
in  load. 

For  the  case  of  a neat  fit,  use  of  the  load  distribution  in  our  finite  element  model  has  given 
results  that  compare  reasonably  with  experimental  values.  Tangential  strains  are  within  7“„ 
of  the  experimental  values;  however  there  are  larger  discrepancies  in  the  radial  strains.  Stress 
concentration  factors  obtained  from  photo-elasticity,  reference  (4),  are  in  general  agreement 
for  cases  of  clearance,  neat  and  interference  fits.  Specific  values  of  stress  concentration  factors 
for  clearance  and  interference  fits  for  the  geometry  of  the  present  lug  are  not  given;  however 
the  trend  is  toward  higher  stress  concentrations  with  increase  of  clearance  and  lower  stress 
concentrations  for  interference  fits.  For  the  case  of  a neat  fit  the  experimental  value  of  stress 
concentration  is  within  10",,  of  the  theoretical  value 

The  effect  of  various  fits  on  the  stress  distribution  is  now  considered  for  a constant  pin 
load  of  50,000  N.  Consider  firstly  the  stress  distribution  shown  in  Figure  1 1 for  the  case  of  a 
neat  fit.  The  maximum  stresses  are  the  tangential  stresses  on  the  inside  of  the  lug  hole  at  boundary 
node  9 (tl  = 90  ).  Comparison  with  Figure  12  for  the  case  of  a clearance  fit  shows  that  the 
tangential  stresses  at  90  are  greater  for  the  clearance  than  for  the  neat  fit  while  the  tangential 
stresses  at  0°  and  45°  (are  smaller  for  the  clearance  than  for  the  neat  fit).  Radial  stresses  at 
l>  = 45  and  = 0 are  similar  for  the  clearance  and  neat  fit  cases,  though  not.  of  course,  at 
fi  = 90  . 


Comparing  now  the  stresses  for  the  neat  fit  (Fig,  II)  with  the  stresses  resulting  from  an 
interference  of  —0  04  mm  (Fig.  13),  it  is  apparent  that  the  tangential  stresses  at  0 = 90  with 
the  interference  fit  are  lower  than  those  with  the  neat  fit,  while  the  tangential  stresses  at  0 and 
45°  are  slightly  greater  with  the  interference  fit  as  were  the  radial  stresses.  Over  all.  the  effect 
of  —0  04  mm  interference  is  to  reduce  the  maximum  stresses.  However  for  higher  interference 
fits,  as  shown  in  Figures  14  to  16,  progressively  all  stresses  are  increased  including  the  tangential 
stresses  at  90°. 

From  the  point  of  view  of  minimising  the  stresses  in  the  lug  for  a given  pin  load  it  seems 
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that  the  smaller  interference  fits  are  optimum.  A plot  of  maximum  tangential  stress  versus 
interference  is  shown  on  Figure  19.  However,  in  the  case  of  fatigue  it  is  the  increment  of  stress 
during  a stress  cycle  that  is  important.  The  increment  in  tangential  stress  Aor  has  been  calculated 
in  Appendix  IV  and  also  plotted  in  Figure  19,  this  corresponds  to  a stress  cycle  starting  from 
zero  pin  1 /ad.  As  can  be  seen  from  Figure  19  the  minimum  value  of  this  increment  corresponds 
to  higher  interference  fits  of  —0  07  mm;  for  fatigue  it  seems  that  the  higher  interference  fits 
are  optimum. 

Factors  determining  the  accuracy  -'f  the  theoretical  load  distribution  are  now  considered. 
In  part,  the  accuracy  is  dependent  on  the  number  of  points  used  in  determining  the  load  dis- 
tribution. For  clearance  fits  the  small  contact  arc  may  result  in  only  7 or  8 points  determining 
the  shape  of  the  load  distribution.  Also  the  arc  length  between  adjacent  points  is  relatively 
large  in  comparison  with  the  contact  arc.  This  makes  convergence  of  the  solution  more  difficult 
since  jumps  in  contact  length  from  one  point  to  the  next  have  a substantial  effect  on  the  shape 
of  the  load  distribution.  For  neat  and  interference  fits  the  longer  contact  arc  ensures  that  more 
points  are  used  to  determine  the  solution  and  this  improves  the  accuracy. 

The  accuracy  of  the  load  distribution  is  also  determined  by  the  accuracy  of  the  individual 
terms  /tij  comprising  the  matrix  of  influence  coefficients  and  subsequent  inversion.  In  this  analysis 
the  values  of  these  terms  were  only  taken  to  four  significant  figures.  This  was  just  acceptable; 
however  six  significant  figures  is  suggested.  Fortunately,  no  difficulty  was  encountered  in  achieving 
an  accurate  inversion  of  the  matrix  of  influence  coefficients.  In  this  analysis  solutions  would 
not  converge  for  levels  of  interference  above  —0-07  mm  for  loads  of  50,000  N or  for  loads 
lower  than  30,000  N with  small  amounts  of  interference. 

Loss  of  accuracy  may  also  occur  in  the  program  in  equations  (6)  and  (13)  where  the  result 
is  the  difference  between  two  almost  equal  numbers.  In  this  program  satisfactory  results  were 
obtained  using  single  precision  arithmetic.  The  program  was  run  a PDP-10  computer  where 
the  word  length  is  36  Bits  and  gives  at  least  eight  significant  figures  for  single  precision.  For 
computers  with  wordlengths  less  than  this,  double  precision  is  adv^safi'c. 

6.  CONCLUSION 

A method  has  been  developed  for  the  solution  of  the  non-linear  pin-lug  contact  problem. 
This  has  allowed  a finite  element  analysis  of  a pin  loaded  lug  which  has  provided  the  stress 
distribution  around  the  lug.  The  method  is,  however,  restricted  to  the  concept  of  a flexible  lug 
loaded  by  a rigid  pin.  The  example  used  of  an  aluminium  alloy  lug  loaded  by  a steel  pin  appears 
to  fit  within  this  concept,  as  a comparison  between  experimental  and  theoretical  result  shows 
agrement  within  7-10%. 
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APPENDIX  I 


The  data  file  INPUT  requires  the  following; 
Variables 


NPOINT,  NLL,  NITER 
RPIN,  F,  ANGLE,  QS,  THICK 
PLL(\) 

PS{I) 

XA  (/,  J) 

YA  (/,  J) 


Format 


5F12-5 


6EI2-5 
6EI2  5 


NPOINT 


NITER 


ANGLE 


THICK 

PLL(\) 

PS(I) 

XA  (/,  J) 
YA{I,J) 


is  the  number  of  points  around  half  the  lug 

Set  to  1 ; this  was  left  in  to  allow  several  load  cases  to  be  considered  at 
once;  however  the  er.lra  programming  has  not  yet  been  implemented. 

is  the  number  of  iterations  to  be  made 

is  the  radius  of  the  pin 

is  the  fit  of  the  pin  in  the  lug.  negative  for  interference 
is  the  angle  in  degrees  between  adjacent  points 
is  the  total  pin  load  corresponding  to  the  sine  load  distribution, 
is  the  thickness  of  the  lug 

Set  this  to  the  required  pin  load  for  the  complete  lug. 

Radial  force  at  point  I due  to  sinusoidal  loading 
X displacement  at  point  / due  to  a unit  radial  load  at  point  J. 

Y displacement  at  point  I due  to  a unit  radial  load  at  point  J. 


APPESDIX  ii- 


BtSrAVAllABlE  COPX  ~ 


C OROGHAM  TO  SOLVE  CONTACT  PROBLEM  BE  I WEEN  PIN  AND  LUC 

C 

C 

C0MM0N/BL0CKl/NPnlNT,NLL.NlTER.KKK,PLL(13) .scale 

COMMON/PR  I NT R/ JSC A, JSCB.JSCC 

integer  TAPNM1(2) .TA0NM2(2) ,TAPNM3(2) 

data  TAPNM1/6HINPUTS/, TAPNM2/6HCHECKS/,TAPNM3/7H0UTPUTS/ 
data  JSCA/1/, JSCB/2/, JSCC/3/ 

0PEN(UNIT  = JSCA,FILE«TAPNM1, ACCESS* 'SEUIN'  ) 

OPEN < UN  I T = JSCB.nLE«TAPNM2,  ACCESS* 'SEQOUT'  ) 
0PEN<UNIT=JSCC.riLE»TAPNM3, ACCESSs’SEUOUT ' ) 

0PEiM(UNIT  = 4,riLEs’0I3OAT$'  ) 

scale*, la 
C 

CALL  read 

c 

CALL  SINE(NPOINT) 

C 

CALL  COORU 
C 

CALL  MATA 
C 
C 

DO  laac  KKK=1, NITER 


C 

500 

C 

C 

C 

C 

C 

C 

1000 


10 

20 

30 

C 

C 

c 


c 

c 

c 

c 


CALL  PniNC 

call  rad 

CALL  RLOAU 

CALL  TLOAU(I) 

IF( I ,EQ. 1)GOTO  500 

CALL  DISP 

CALL  WRITE 

CONTINUE 

END 

SUBROUTINE  READ 

COMMON/BLOCKI/NPOINT.NLL, NITER.KKK,PLL(10) .SCALE 
COMMON/BlUCK2/XA(20,20) ,YA(20,20).A(20.20) 
COMMON/BLOCK4/XOlSPS(2a) , YD  I SPS ( 20 > » PS ( 20 ) 
COMMUN/BLUCKA/RPIN.F, angle. QS.QT.R( 20) .THICK 
COMMON/PR  1 IN  TR/ JSC  A,  JSC  3.  JSCC 
F0RMAT(3I5) 

F0RMAT(7F12.5) 

F0RMAT(6E12.5) 

READ  IN  number  of  POINTS,  ’IN  LOAD  AND  NUMBER  OF  ITERATIONS 

REAO< JSCA.10)NPOINT, nLL'NITER 
WRITE! JSC«, 10) NPOl NT, NLL, niter 

READ  IN  RADIUS  OF  PIN,  FIT  ,PIN  LOAU,  OS  AND  THICKNESS  OF  LUC 
FIT  IS  NEGATIVE  FOR  INTERFERENCE, 


READ! JSCA,20)R?r),F, angle, OS, thick 


BEST  AVAILABIF  COPY 


c 

c 

c 


c 

c 

c 

c 

c 

100 

c 

c 

c 


200 


c 

c 

c 

c 


100 


200 


C 

C 

c 


c 

c 

c 


100 


^RITtlC  JSCH,23)RPIN,r,  asigle:»qs»thick 

READ  IN  particular  pri»  LOADS 

REAOt JSCA»30) (PLLI I > » I*1#NLL) 

WRITE!  JSCb,.3a)  (PLL(  I > . I-l.NLL) 

DO  100  IbI.NPOINT 

READ  IN  radial  LOADS  FOR  SINUSOIDAL  LOADING  3S 
< OR  ANY  other  approximate  LOADING  I.E,  CONCENTRATED 
LOAD  AT  VERTEX  OF  LUC  WILL  WORK  ) 

READ! JSCAi30)PS(  I > 

WRITE! JSCB.30)PS!  I ) 

READ  IN  X AND  Y DISPLACEMENTS  OF  EACH  POINT  DUE  TO  UNIT  LOADS 

DO  200  Isl»NPOINT 

READ! JSCA.30) !XA! I ,J) , jal,NPOINT> 

WRITE! JSCB,30) (XA! 1 , J) , J«1,NP01NT) 

READ! JSCA.30) ! YA! I , J) , J«l»NPOINT) 

WRITE! JSCH,30) !YA! I , J) , JsI.NPOINT) 

RETURN 

END 

subroutine  SINE!NP0INT) 

COMMUN/BLOCK2/XA!20,20) ,YA!20»20) . A!20,20) 
COMMUN/BLOCK4/XDISPS!20) , YDISPS!20> •PS!2a) 

FORM  displacements  CORRESPONDING  TO  RADIAL  SINE  LOAD 
OR  any  other  APPROXIMATE  LOAD  DISTRIBUTION 

DO  200  Ial»NPOINT 
SUMX=0. 

SUMY»0, 

DO  100  Jsl.NPOINT 
SUMXaSlIMX^XA!  I , J>*PS!  J) 

SUMY  = S1)HY*YA!  I » J)*PS!  J) 

XOISPS! I )=SUMX 

YDISPS! I )=SUMY 

CONTINUE 

RETURN 

END 

subroutine  COORD 

COMMUN/BlOCK1/NPOINT,NLL.NITER»KKK,PLL! 10) .scale 
COMMUN/BLOCK3/XCOORD!20) , YCOORD!20) 
C0MM0N/BL0CK5/RPIN,F,ANGLE,QS.QT.R! 20) .THICK 

RLUG  IS  RADIUS  OF  HOLE 

rlug=rpin*f 

THETAS0. 

NOW  CALCULATE  COOROS  OF  UNDEFLECTED  LUG 

DO  100  Isl.NPOINT 
XCOORD!I)=  RLUG*SIND! T NETA) 

YCUURO! I ) =-RLUG*COSO! THETA) 

THETAsTHEU^ANGLE 

RETURN 

END 


r ’ — m\  mum  iufy — ^ 

SUBROUTInL  MATA 

COMMON/BLOCKl/NPOlNT,NLL»NITER.KKK»PLL(li’)  »SCAlE 
COMMON/0LOCK2/XA(2IZ),20)  , YA(20»20>  »A(2K*20) 

COMMON/BLOCK3/XCOORD(2a) . YCOORO(20> 

C 

c matrix  a is  a set  or  radial  displacements  for  unit 

c radial  loads, 

, c 

DO  200  lal.NPOINT 
* 00  200  J»1.NP0INT 

c 

c WE  MUST  NOW  DETERMINE  WHICH  RADIAL  DISPLACEMENTS 

C ARE  POSITIVE  AND  NEGATIVE 

C radial  compression  is  taken  to  be  POSITIVE, 

c R IS  radial  vector  to  point 

C RDIS  IS  RADIAL  VECTOR  TO  DISPLACED  POINT 

C 

R»SuKT( (XCOOROC I > ) •*2^ ( YCOORD ( I ) )**2) 

RDIS*SQRT( (XCOOROt I )*XA( I , J) YCOOROI I J*YA( I , J) )»*2) 

A( I , J)»ROIS-R 
200  CONTINUE 

CALL  GENINV{A»NPOINT,20.O» IRRORfCONU) 

return 

END 

SUBROUTINE  GEN  I NV ( A , N , M , D » I RROR * CONO ) 

C 

C this  INVERSION  ROUTINE  WAS  WRITTEN  bY  D,W,G,  MOORE 

C COMPUTING  CENTRE,  UNIVERSITY  OF  WESTERN  AUSTRALIA 

i C 

DIMENSION  AIM.M) , IPIV(100) , INO(100»2> 

IF(N,LE .0.OR.N,GT.M,OR,M,GT.100)  goto  80 
CALL  OVERTL(JJJ) 

ASSIGN  7777  TO  G08ACK 
GOTO  7779 

I 7777  0=1,0 

IRROR=n 
DO  10  I=liN 
10  IPIV(I)=0 

DO  220  I=1,N 
AMAX=a,0 

c 

c SEARCH  SUB-MATRIX  FOR  LARGEST  ELEMENT  AS  PIVQT 

C 

DO  70  J=1,N 
IF(IPIV{ J) )102»30»70 
30  DO  60  Wrl.N 

IF( IPIV(K)-l)40,60,ia3 
C 

C THIS  ROW  HAS  been  A ^IVOT 

C 

40  ir( AHS( A( J,K) ) ,LE. AMAXIGOTO  60 

50  IROW=J 

IC0L=K 

AMAX  = ABS( A{ J,K  ) > 

60  CONTINUE 

i 70  CONTINUE 

C 

C PIVOT  FOUND 

C 

1PIV(  IC0L)  = 1PI  V{ ICOLI^I 


BEST  AVAIUBIE  COPY 


90  ir( IHOW.EQ. ICOL)COTO  130 

C 

C 'lAKE  PIVOT  A DIAGONAL  ELEMENT  BY  ROW  INTERCHANGE 

C 

95  Os-D 

DO  100  K*1>N 
AMAX«A( IROW,K) 

A( IR0W,K)»A( ICOL»K) 

100  A( ICDL»K)aAMAX 

130  INO(  I .DalROW 

IND( I ,2)«IC0L 
AMAX»A( ICOL# ICOL) 

IF(D.LT.0.1E-15)O«0, 

IF( AMAX.lT.0,1E-15)AMAXb0, 

D«D*AMAX 

A{ ICJL» ICOL)»1,0 
C 

C DIVIDE  ROW  BY  PiVOT 

C 

DO  1A0  K3l,N 

140  A( lCOLiK)=A( IC0L,K)/AMAX 

DO  220  Jsl.N 
IF( J.EQ. ICOL)GOTO  220 
180  AMAXsA( J, ICOL) 

A< J, ICOL)'0. 

DO  1V0  K«1»N 

190  A( J,K)cA( J,K)-A< IC0L.K)*AMAX 

220  CONTINUE 

C 

C FOR  INVERSE  OF  A,  INTERCHANGE  COLUMNS 

C 

230  DO  260  Isl.N 

JsN+l-I 

240  IROWsIMDE J.l) 

IF< INO( J,l) .EQ. IN0( J,2) )COTO  260 
IC0L=IND( J,2) 

DO  220  Ksl.N 
AMAXaAEK. IROW) 

A{K, IROW)*A(K» ICOL) 

250  A(K,  ICODsAMAX 

260  CONTINUE 

270  CALL  OVERFL(JJJ) 

IF( JJJ,NE.2) IRROR'l 
ASSIGN  6000  TO  GOBACK 

consav«cono 

GOTO  7779 

8000  CONO=CONSAV*CONO 

RETURN 

80  IRR0K=2 

return 

102  IRR0Ra3 

return 

103  CONTINUE 
IRHORsA 

return 

7779  CONO»0, 

DO  7/80  1*1, N 
DO  7780  Jsi.N 

7780  CONOaCONO^A( I , J)»*2 
CONO=SORT(CONO) 


L ^ J 


iiOTO  GOBACK 
ENJ 

subroutine.  POINC 
COMMaN/BLOCKl/NPOINT,NLL»NITER»KKK,HLL(Xa) » scale 
COMMJN/BLOCK4/XOISPS(20)  » YDlSPS(2k))  *PS(2a) 
COMMr)N/BLCJCK5/DISPX(20>  ,DISPY(20)  ,P(20) 

THIS  ROUTINE  updates  LOADS  AND  DISPLACEMENTS 

IF(KKK,GT.1)G0TD  200 
DO  100  I=l.NPOINT 
OISPX( I )5SCALE*XDISPS(  I ) 

0ISPY{ I )sSCALE»YdISPS( I ) 

P( I )»SCALE*PS(  I > 
return 

END 

subroutine  RAD 

COMMnN/BLOCKl/NPOlNT,NLLiNITER.KKK»PLL(10) iScALE 
C0HM0N/BL0CK3/XC00RD(2a) ,YCOORD(20) 

COMMON/BLUCK5/OISPX(20) . D I SPY ( 20 ) , P ( 20 ) 

C0HM0N/BL0CK6/RPIN,F, angle, QS»QT.R( 20 > .THICK 
IF(F,GT.0. )G0T0  140 
IFCF.LT .3. )G0T0  350 

NEAT  FIT  - contact  LENGTH  OVER  90  DEGKEES 

ICL=10 
GOTO  550 

FIND  CONTACT  LENGTH  FOR  CLEARANCE 

1 = *) 

I = I>1 

IF(P( I ) .EQ.0. )G0T0  250 

IF( I ,EQ.NP0INT)G0T0  200 

GOTO  150 

icl=npoint 

GOTO  330 

ICL=1*1 

SHOULD  BE  ICL=l  . HOhEVER  ICL»I*1  ALLOWS  SOLUTION  TO  CONVERGE 
NO  error  in  THIS  ASSUMPTION  PROVIDED  SOLOTIUN  DOES  CONVERGE 

CONTINUE 
GOTO  550 

FIND  contact  length  FOR  INTERFERENCE 

l-o 

1-1*1 

IF(P( I ) .EO.0. )GOTO  500 

IF(  I .EQ.NPOI.NT)GOTO  450 

GOTO  400 

ICL=NPniNT 

GOTO  550 

ICLsl^l 

CONTINUE 

FORM  RADIAL  DISPLACEMENTS  THAT  ARE  10  BE  APPLIED  TO  LUG 


BtSfAVAllABlE  COPY 

DRPr'JsKPIN 
00  1130  I=1,NP0INT 
THETAaTHETA4.ll, 25 

RH=SURT{ (XCOORO( I )4DISPX( I ) )**2*(YCOOHD(I )*0ISPY( I ) 
-DISPY(1)4F)**2) 

R(I  )>'RPIH-RH 

REHOVE  tensile  RADUL  DISPLACEMENT  CORRECTIONS 
OUTSIDE  contact  LENGTH 


IF(  I ,LT.  ICDGOTO  100 
IF(RPIN.LT.RH)R( I )»0. 

CONTINUE 

RETURN 

END 

SUBROUTINE  RlOAD 

common/blocki/npoint,nll.niter.kkk,pll(i2) » scale 

COMMON/BLUCK2/XA(20,20) » YA(20»20) ,AC20|20) 
COMMON/BLOCK5/OISPX(20) *DISPY(20) ,PC20) 
COMMON/BLOCK 6/RP I N, F, angle, QS.QT,R( 20) , THICK 


NOR  FORM  RADIAL  LOADS 
ISHa0 

DO  200  I=1,NP0INT 
SUM=0. 

00  100  J=1,NP0INT 
SUM  = SIJM4.A(  I , J)*R(  J) 

REMOVE  LOADS  FOR  WHICH 
radial  displacements  are  iJERO 

IF(R(  I ) .F(J.0,  )SJM  = 0, 

ADO  INITIAL  LOAD  TO  OBTAIN  TOTAL  LOAD 

P(  n=P(  I )4SUM 

REMOVE  tensile  LOADS 

IF(P( I ) .LT.0. )P( I )=0. 

IF(P< I ) .Ey.0.  ) ISW  = 1 
IF(  IiiW,EQ.l)P{  I )=0. 

CONTINUE 

RETURN 

END 

SUBROUTINE  TLOAD(K) 

COMMnN/BLOCKl/NPOINT,NLL,NITER,KKK,PLL(l^) iSCALE 
COMMON/BLOCK5/OISPX(20) ,DISPY(20) ,P(20) 
C0MMUN/BLUCK6/RP1N,F, angle, QS,QT,R( 20) .THICK 

THIS  ROUTINE  ADJUSTS  THE  LOAD 
DISTRIBUTION  TO  GIVE  THE  REQUIRED  TOTAL  LOAD 

K = 0 

ICOUNTaS 

ICOUNTalcuUNT+l 

THETAaC. 

OTapCl) 


■ ^VMIABU  lur  i ^ 

DO  11^0  uiiiNPorn  ‘ j 

THETA  a THETA^ANGLE 
aT3QT4P( I )•( COSD( THETA) ) 

100  CONTINUE 

IF(KKK,GT.1)G0T0  250 
C 

C THE  THEORY  HERE  IS  THAT  WE  ARE  CLOSE  ENOJCHT  TO  THE  LOAD 

C TO  MAKE  AN  ASSJMPTION  OF  LINEAR  LOAD  WITH  LOAD  OISTHIBUTION 

C RELATIONSHIP 

C 

IF(F,LT .0. )GOTO  350 
C 

C 'JSE  A FACTOR  X FOR  NEAT  OR  CLEARANCE  FITS 

C 

250  XcABS{PLL(1)/(0T*2.  ) ) 

DO  200  I=l»NPOINT 
200  P(I)=P(I)*X 

IF( ICOUNT.GT,100)GOTD  300 
IF{X.LT.0.99,OR.X.GT,1,0X)GOTO  50 
300  RETURN 

C 

C JSE  FACTOR  SCALE  FOR  INTERFERENCE  FITS 

C 

350  CONTINUE 

aTT=2 . *QT 

IF(QTT,LT.PLL(1) )G0TD  400 
C 

c take  first  value  greater  Than  plluj 

c 

GOTO  303 
C 

C SCALE  IS  NOW  SET  FOR  ALL  OTHER  ITERATIONS 

C 

400  IF(KKK,GT.1)G0T0  300 

scale=scale*.i 

IF(Sc:ALE.GT.100.  IGOTO  303 
K = 1 

RETIJKN 

END 

SUJRUUTINE  DISP 

COMMON/BLUCKl/NPniNT,NLL,NITER.KKK,HLL(10) .scale 
COMMC)N/aLUCK2/XA(20,2Z)  , YA(20,20)  .A(20,20) 

COMMON/BLUCK5/DISPX ( 20).DISPY(20).Pi20) 

C 

C form  SET  OF  displacements  CORRESPONDING  TO  LOAD  DISTRIBUTION 

C 

DO  200  lal.NPOINT 
SUMX=0, 

SUMY=0, 

DO  100  J=1,NP0INT 
SUMX  = SUMX'»XA(  I . J)*P(  J) 

100  SUMYsSUMY-.'YAI  I , J)*P(  J) 

DIbPX( I )=SUMX 
200  DISPY( I )=SUMY 

return 

END 

subroutine  write 

commjn/blocki/npoint.nll,niter,kkk,rll(13) .scale 

CO'1MON/BLOCK5/DISPX(20)  .0ISPY{23)  .P(20) 

C0HMJN/BlUCK6/RPIN,F , ANGlE,QS.OT,H(20) ,THICK 


BtSI  AVAimBlt  m — ' 


common/PR  luT^/jsc  A,  j'ics.jscc 

OI^E'^SION  AP(2il)  , APX(23)  , APY(20)  ,APKESS(2id)  ,rX(40>  ,FYU0) 
<MI=1000 

IF(KKK.LT.KMI )GOTO  230 

10  FORMATC  iteration  M4,12X, 'FORCE.  N,',12X, 

1 'DISPLACEMENT  MM  12  X ,' CONTACT  PRE.SSURE  MPA,'//) 

20  FORMATC  PIN  LOAD  ' , 2X  , F 10 . 2 , 20X  , * FIT  ' » 2X . F10 . 6 » 20X , 

1 ' RPIN  ' ,2X,F13.6//) 

40  FORMATC  POINT' ,5X, 'RLOAD* ,9X, 'XLOAU  '.BX,'  YLOAO  '.SX, 

1 ' XDISPL'»3X,’  YOISPL' ,4XCC0NTACT  PRESSJRE'//) 

50  format (2X. I4,F13.2,2(2X,F13,2) ,2(2X»E13,6) ,2X,F13.2) 

60  FORMaT(1H  ///) 

C 

C FI  and  F2  are  conversion  FACTORS  FOK  UNITS 

C 

Fls4 ,44622 
F2»25.4 

AP(l)=ri*(2.3«P{l)  ) 

DO  90  I=2*NP0I  JT 
90  AP( I )=F1*(P(  I ) ) 

DTTsF1*(2.0*DT) 

T=F2»(THICK) 

RL  = F2»(RPIN-^F) 

c 

C NOW  FORM  CONTACT  PRESSURE 

C 

3 = 57.296/(RL*ANGlE.*T) 

DO  95  I=l.NPOINT 
95  APRESS( I )=AP( I )*S 

WRITE!  JSCC.iaXKK 
FIT=F2*(F) 

RP=F2*(PPIN) 

WRITE!  JSCC,  22)  :jTT.  FIT,  RP 
WRITE! JSCC, 43) 

THLTA  = e . 

DO  100  I=l,NPOI'JT 
APX(I)=  AP! I )*5InO!ThETA) 

APY!I)=-AP!I)*COSO(ThETA) 

theta=theta*angle 

DX  = F2*!DISPX!  I ) ) 

DY  = F2*!0ISPY(  I ) ) 

100  WRITE! JSCC, 53) I • A^! I ) , APX! I ) , APY! I ) * DX , OY , APRESS C I > 

WRITE! JSCC, 63) 

200  return 
END 


APPENDIX  III 


Stress  concentration  factors 

Referring  to  Figure  4 with  er^  being  the  tangential  stress  at  point  9,  then  the  stress  con- 
centration factor  is  given  by; 


A 


where  PI  A is  the  nominal  stress 

i Referring  to  the  dimensions  given  in  Figure  3,  the  minimum  cross  sectional  area  of  the  lug 

' is  given  by: 

' A = 2 X (29-591  - 16-866)  x 12-700  x 10  " m* 

= 3-232  X 10-‘m* 

For  a pin  load  of  50,000  N,  the  nominal  stress  is  PjA  = 1 -547  x 10*  Pa. 

For  the  pin  load  of  50,000  N the  following  stress  concentration  factors  are  given  for 
various  fits. 


Fit  (mm) 

a*  (MPa) 

K 

0-05 

395-9 

2-559 

O-OO(Neat) 

333-4 

2-155 

-0-04 

230-8 

1-492 

-0-05 

236-7 

1-530 

-0-06 

252-6 

1-632 

-0-07 

279-4 

1-806 

The  experimental  value  of  K for  a neat  fit  given  in  reference  (4)  is  approximately  2 -4. 


APPENDIX  IV 


Tangential  Stresses  due  to  Interference  only,  zero  pin  load 

From  reference  (6)  the  tangential  stress  at  the  pin-lug  interface  due  to  an  interference  pin 
with  zero  pin  load  is  given  by: 


and 


where 


OT 

Ee 

9 

Ee 


Eu 

E, 

d 


H’ 


Vi  are  Youngs  modulus  and  Poisson’s  ratio  for  pin 
V are  Youngs  modulus  and  Poisson’s  ratio  for  lug 
is  diameter  of  lug  hole 
is  width  of  lug 


e = 2Fld. 

Here  we  take  E^  = 206-8  x 10®  MPa  and  vj  = 0-30,  while  E and  v have  the  same  values  as 
previously  given. 

Substituting 

ar  = 32-75  X 10®  X F MPa 


Interference  (mm) 

OT  (MPa) 

-0-04 

131-0 

-0-05 

163-7 

-0-06 

196-7 

-0-07 

229-3 

The  following  maximum  values  of  ot  have  been  extracted  from  tables  13-20;  from  these 
values  Aar  is  calculated.  Aar  is  given  by: 


Aar  = <^7’max  -~  <rr  (zero  pin  load) 


Fit  (mm) 

<rrmax(MPa) 

Aar  (MPa) 

0-05 

395-9 

395-9 

0-00 

333-4 

333-4 

-0-04 

230-8 

99-8 

-0-05 

236-7 

73-0 

-0-06 

263-3 

66-8 

-0-07 

294-9 

65-6 
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TABLE  13 


Stresses  in  Lug  with  Neat  Fit  Pin ; Pin  Load  of  50,000  N 


Radial 

Distance 

(mm) 

Radial  line  through 
Point  1 

Radial  line  through 
Point  5 

Radial  line  through 
Point  9 

OT 

OHT 

<^R 

O'f 

ORT 

OR 

O'P 

ORT 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

16  86 

-107-2 

rm 

-3-8 

-125-1 

159-6 

-1-2 

-42-8 

333-4 

- 4-8 

17-62 

- 94-3 

IBKI 

-0-7 

-110-1 

157-9 

-3-5 

-34-6 

294-2 

- 8-0 

18-38 

- 86-0 

142-9 

-1-3 

-101-0 

154-3 

-5-8 

-22-0 

257-8 

-16  0 

19-65 

- 70-3 

141-5 

-0-9 

- 80-4 

151-8 

-8-3 

-13-6 

218-7 

-16-2 

20-92 

- 59-0 

139-8 

-1-1 

- 67-0 

146-3 

-9-8 

- 5-2 

183-5 

-17-4 

22-70 

- 42-8 

141-4 

-1-0 

- 46  -5 

144-0 

-9-6 

- 4-3 

149-9 

-12-8 

24-48 

- 30-3 

142-4 

-1-0 

- 31-6 

139-6 

-8-8 

- 0 8 

121-4 

-10-3 

27-05 

- 13-9 

146-5 

-0-5 

- 13-9 

140-1 

-4-9 

- 1-4 

84-2 

- 4-7 

29-59 

0-5 

151-4 

-0-7 

0-0 

142-6 

-1-5 

1-4 

48-3 

- 0-8 

TABLE  14 


Stresses  in  Lug  with  0.05  mm  Clearance  Fit  Pin ; Pin  Load  of  50,000  N 


Radial 

Distance 

(mm) 

Radial  line  through 
Point  1 

Radial  line  through 
Point  5 

Radial  line  through 
Point  9 

OR 

(77* 

ORT 

OR 

OT 

ORT 

OR 

OT 

ORT 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

16-86 

-118-1 

-1-8 

-133-9 

137-7 

-9  0 

395-9 

-4-7 

17-62 

-103-3 

106-7 

-1-7 

-117-1 

140-9 

-0-7 

-1-4 

18-38 

- 96  5 

111-9 

-1-6 

-109-0 

140-9 

13-8 

289-7 

-7-4 

19-65 

- 78-7 

120-4 

-1-6 

- 86-6 

143-0 

-13-0 

14-4 

235-5 

-5-2 

20-92 

- 66  -8 

126-5 

-1-5 

- 73-0 

140-3 

-15-8 

19-5 

190-0 

-7-5 

22-70 

- 48-2 

135-8 

-1-4 

- 49-9 

140-5 

-16-0 

116 

146-8 

-3-7 

24-48 

- 34-2 

143-8 

-1-2 

- 33-7 

137-6 

-14-8 

9-7 

110-3 

-3-1 

27-05 

HSI 

156-6 

-0-7 

- 14-6 

140-8 

- 8-1 

2-6 

60-6 

-1-2 

29-59 

■ 

170  0 

-0-7 

0-1 

146-6 

- 1-8 

1-4 

10-7 

-0-8 

TABLE  15 

Stresses  in  Lug  with  —0.04  mm  Interference  Fit  Pin;  Pin  Load  of  50,000  N 


Radial 

Distance 

(mm) 


Radial  line  through 
Point  I 


Radial  line  through 
Point  5 


OT  ORT 


Radial  line  through 
Point  9 


-114-4 
-95-9 
-85-9 
-65-0 
-51-1 
-33-8 
-20-4 
- 8 6 
0-8 


TABLE  16 

Stresses  in  Lug  with  -0.05  mm  Interference  Fit  Pin;  Pin  Load  of  50,000  N 


Radial 

Distance 

(mm) 


Radial  line  through 
Point  I 


OR  OT 

(MPa)  (MPa) 

-134-8  232-6 

-115-2  221-3 

-103-6  208-1 

- 82-0  195-4 

- 66-9  181-9 

- 46-8  171-1 

- 31-4  159-8 

- 14-1  151-4 

0-5  145-6 


Radial  line  through 

Radial  line  through 

Point  5 

Point  9 

OR  OT  ORT  OR 

(MPa)  (MPa)  (MPa)  (MPa)  (MPa)  (MPa) 

-132-5  236-9  -1-8  -116-3  236-7  -16 

-114-7  225-7  -1-1  -99-0  224-7  0-0 

-103-0  213-1  -2-1  -88-1  212-0  -0-4 

- 80-7  199-1  -2-2  - 67-5  197-4  0 5 

- 65-5  184-2  -2-7  -53-3  182-2  0-7 

- 45-4  171-7  -2-3  -35-8  169-1  2-5 

- 30-7  158-9  -2-4  -21-7  155-8  3-0 

- 13-5  148-6  -1-4  - 9-2  147-9  5-0 

0-0  139-8  -1-4  0-7  144-7  5-2 


TABLE  17 


Stresses  in  Lug  with  —0.06  mm  Interference  Fit  Pin;  Pin  Load  of  50,000  N 


Radial 

Distance 

(mm) 

Radial  line  through 
Point  1 

Radial  line  through 
Point  5 

Radial  line  through 
Point  9 

G'T 

ORT 

OR 

Gt 

Out 

OR 

Gf 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

16-86 

-141-8 

263-3 

-143-1 

260-0 

-2-1 

-121-9 

252-6 

-1-7 

17-62 

-121-4 

248-1 

-123-9 

247-2 

-1-1 

-104-2 

239-  . 

1-0 

18-38 

-108-5 

231-5 

-111-1 

233-0 

-2-2 

-92-1 

225-3 

1-2 

19-65 

- 85-8 

215-3 

-1-3 

- 87-0 

216-9 

-2-0 

-71-0 

209-7 

3-1 

20-92 

- 69-7 

198-8 

-1-9 

-70-5 

200-2 

-2-5 

-55-8 

194-1 

3-9 

22-70 

- 48-9 

185-3 

-1-2 

- 48-8 

185-8 

-2-1 

-37-8 

181-2 

6 0 

24-48 

- 32-8 

171-5 

-1-5 

- 32-9 

171-3 

-2-3 

-22-9 

168-5 

6-5 

27-05 

- 14-7 

160-3 

- 14-5 

159-3 

-1-4 

- 9-7 

162-1 

7-5 

29-59 

0-6 

1 

152-0 

H 

0-0 

148-9 

-1-5 

160-9 

6-4 

TABLE  18 


Stresses  in  Lug  with  — 0.07  mm  Interference  Fit  Pin ; Pin  Load  of  50,000  N 


Radial 

Distance 

(mm) 

Radial  line  through 
Point  1 

Radial  line  through 
Point  5 

Radial  line  through 
Point  9 

<r/i 

GT 

ORT 

OT 

G'p 

ORT 

OR 

GT 

ORT 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

(MPa) 

(MPa)  ' 

1 

(MPa) 

(MPa) 

16-86 

-156-0 

294-9 

-6-8 

-158-1 

290-4 

-2-3 

-132-7 

'>79-4 

-1-9 

17-62 

-133-6 

217-2 

-1-4 

-136-8 

275-7 

-1-2 

-113-6 

264-2 

1-5 

18-38 

-119-2 

258-2 

-2-6 

-122-5 

259-5 

-2-4 

-100-1 

248-5 

2-0 

19-65 

- 94-2 

239-5 

-1-4 

- 95-9 

241-2 

-2-2 

-77-2 

230-9 

4-6 

20-92 

- 76-5 

220-5 

-2-1 

- 77-5 

222-0 

-2-8 

-60-5 

213-6 

5-8 

22-70 

- 53-5 

204-9 

-1-3 

- 53-7 

205-6 

-2-3 

-41-0 

199  6 

8-3 

24-48 

- 35-9 

189-0 

-1-7 

- 36-1 

188-9 

-2-5 

-24-7 

185-8 

8-9 

27-05 

- I6-I 

175-8 

Bll 

- 15-9 

174-9 

-1-5 

-10-5 

179-5 

9 4 

29-59 

0-6 

166-0 

m 

0-0 

162-8 

-1-7 

0-8 

178-9 

7-4 

TABLE  19 

Experimental  Strains  in  Lug  for  Neat  Fit  Pin;  Pin  Load  of  SO, 


Gauge 

Number 


Gauge  Location 

Type  of 

Value 

Distance  from 
centre  of  Lug 
hole  (mm) 

Angular 

coordinates 

degrees 

strain 

measured 

of  strain 
10'“ m/m 

18-14 

90 

Tangential 

4-12 

18-42 

45 

2-25 

18-24 

0 

2-25 

28-38 

90 

0-67 

27-95 

45 

*» 

2-14 

27-88 

0 

1-86 

17-63 

90 

Radial 

-2-48 

17-77 

45 

-1-80 

17-64 

0 

-1-20 

TABLE  20 

Theoretical  strains  in  lug  for  neat  Fit  Pin;  Pin  load  of  50,000  N 


Radial 

distance 

(mm) 


Radial  line 
through  point  1 


fR 

10“’ m/m 


€T 

10“  ® m/rn 


Radial  line 
through  point  5 


-2-481 
-2-262 
-2-117 
-1-817 
-1  604 
-1-304 
-1-075 
-0-827 
-0-642 


Radial  line 
through  point  9 


-XCOORDj 


A 


FIG.  5 RADIAL  DISPLACEMENT  OF  POINT  i DUE  TO  UNIT  RADIAL  LOAD  AT  j 


Deflected  shape  of  lug  hole  relative  to 
the  pin  (After  vertical  translation) 


FIG.  6 GEOMETRY  OF  PIN  AND  LUG 


FIG.  9 LOAD  DISTRIBUTIONS  FOR  A CLEARANCE  FIT  OF  0.05mm 


FIG.  10  LOAD  DISTRIBUTIONS  AROUND  LUG  FOR  AN  INTERFERENCE  FIT  OF  - 0.04mm 


_2ool-  Radial  line  through  node  1, 0 = 0* 
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FIG.  1 1 THEORETICAL  STRESSES  FOR  A NEAT  FIT  ; PIN  LOAD  IS  50kN 


FIG.  12  THEORETICAL  STRESSES  FOR  A CLEARANCE  FIT  OF  0.05mm  ; 
PIN  LOAD  ISBOkN 
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Stress  MPa 


Radial  line  through  node  5,d  = 45* 


Stress  MPa 
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FIG.  13  THEORETICAL  STRESSES  FOR  AN  INTERFERENCE  OF  -0.04mm  ; 
PIN  LOAD  IS50kN 
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Radial  line  through  node  9,  0 = 90* 


Stress  MPa 


Stress  MPa 
400 1 


300 

200 

100 

0 

-100 


-200 
Stress  MPa 
3001- 
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Stress  MPa 
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200  - 

100- 

0 
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-2001- 


Radial  line  through  node  9,  6 = 90* 


Radial  line  through  node  5,  0 = 45° 


Radial  line  through  node  1,0=0* 
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FIG.  15  THEORETICAL  STRESSES  FOR  AN  INTERFERENCE  OF  - 0.06mm  ; 
PIN  LOAD  IS  50kN 
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Radial  line  through  node  5,  0 = 45 


Radial  lina  through  node  1,0  = 0* 
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FIG.  16  THEORETICAL  STRESSES  FOR  AN  INTERFERENCE  OF  - 0.07mm  ; 
PIN  LOAD  IS50kN 


Radial  lines  along  which 
stresses  and  strains  are 
considered 


FIG.  17  POSITION  OF  STRAIN  GAUGES 


Radial  line  through  node  9,  0 = 90' 


Strain  XIO^ 
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FIG.  18  THEORETICAL  AND  EXPERIMENTAL  STRAINS  FOR  A NEAT  FIT 
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